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Abstract — In this technical note, we present the formulae 
of the derivatives of the Mean Value Coordinates |2| based 
transformations using an enclosing triangle mesh, acting as a 
cage for the deformation of an interior object. 



I. Background 
A. Mean Value Coordinates for Closed Triangular Meshes 

Mean Value Coordinates for closed triangular meshes were 
introduced in |2|. In this section, we briefly review this work 
and describe the notations we will use in the rest of this note. 

As written in 1 2 1, a 3D point rj can be expressed as a linear 
sum of the 3D positions pi of the vertices of a triangular mesh 

For a point x onto the surface (a two-dimensional param- 
eter), we note as usual (t)i[x] the linear function on M that 
takes value 1 on vertex i and on other vertices, and p[x] its 
3D position. The definition of the weights Xi should guarantee 
linear precision (i.e. rj — Ai(ry)pi). 

Ib^{m) bd-rSl ^'^'X^) = (*e integral of the unit 
outward normal onto the unit sphere is 0), we have 



(1) 



Bjj{M) being the projection of the manifold AI onto the unit 
sphere centered in rj. 

Writing that Vx p[x] = J2i 4'i[x]Pi, with J2i — 1, we 
obtain 



The weights Aj are given by 

^ ^ JB^(M) \pIx]-v\"'^v\^) 

Ib„{m) \plx]-,j\'^^n(^) 
And the weights Wi such that ~ ^"''^ are given by 

0j [a; 
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„(M) -v\ 



dSri{x) 
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This definition guarantees linear precision; it gives a linear 
interpolation of the function onto the triangles of the cage; 
and it extends it in a regular way to the entire 3D space. 



Computing the weights Wi: The support of the function 
(t)i[x\ is only composed of the adjacent triangles to the vertex 
i. Then, we can rewrite Eq. |4jas Wi = X^TeiviCi) ^ ^^'^ 

<Pi [x] 



iB^iT) \P[X] -V\ 

Given a triangle T with vertices ^1,^2,^3, we see that 
V wt, ■ [Pt, -V)= I r 1 , dT 

-^B^T) \pm~v\ 



(5) 



„(T) \P[X] - V\ 



(6) 



dT = 



This last integral is simply the integral of the unit outward 
normal on the spherical triangle T. 

By noting nf = with = {pt^^, - 77) A (pt,+2 " v) 

(see Fig|T]l, it can be easily expressed as 



El' 



(7) 



This comes from the fact that the integral of the unit outward 
normal on a closed surface is always 0. 
Finally, we obtain 



(8) 



This point was discussed in p|. As the authors pointed out, 
by noting the 3 by 3 matrix {pt^ — r],pt2 — V^Pts ^ v}^ 
we can derive the weights wf. by 

{wl,wl,wlY =A^-' -m^ (9) 
Since A^^^* • {ptj — r;) = Vi 7^ j, we have from Eq. [sjthat 



■ [Pu ~ V) 



V77 ^ Support (T) 



(10) 




Fig. 1. Triangle T projected on the spherical triangle T. 
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II. MVC Derivatives 

We now present the derivatives of the Mean Value 
Coordinates. Deforming the cage mesh with f{pi) = Pi 
induces a deformation of the 3D space by / = ^ ■ Ai • pi- 
In the rest of the document, for any function h : E ^ F, 
we note dxh,dyh,dzh its derivative by x,y, and z, sjh its 
gradient, Jh its jacobian, and Hh its hessian. 

The deformation function / as defined acts now on 
entirely. The derivatives of / can be expressed as a linear 
sum of positions pi = {xi,yi,Zi}: 



(11) 



Hifv) =T.,m-HK 

Consequently, it allows to specify implicit equations on the 
cage in a linear system by giving specified rotations and scales 
on 3D locations, or to minimize the norm of the hessian to 
force rigidity, as done in the case of Green Coordinates in 
Since A, 



We also have \fc = x, y, 
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dc{wi) ■ J2j V^«i 



2w, 



(E, vw'j 
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or 
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'Mi Y,j HWj 



vw^ ■ E, vK")* + E, vK) ■ V 



(14) 



From these expressions, we see that, in order to get S7Xi{r]) 
and HXi{rj), we first need to obtain \/Wi{r]) and Hwi{ri) for 
each vertex i of the cage. 

Special case: rj lies on the surface of the cage 

Mean Value Coordinates define an interpolation process. 
The function represented onto the vertices of the cage (in our 
case, a space transformation) is extended to the interior of the 
triangles with linear interpolation on each triangle. Then it is 
extended to the space by means of a surfacic integration of 
the function (see Eq. [TJ. 

Since we represent the cage as triangle mesh in the 3D 
case, the deformation function cannot be anytliing more 



than continuous onto the edges of the cage in 3D. Therefore 
Jacobians and Hessians of the deformation cannot be evaluated 
everywhere on the surface of the cage, and we do not provide 
any formula for Jacobians and Hessians of the deformation 
onto the surface of the cage. 



A. Expression of the Jacobians 

In the general case where det(A^) = (pe, - vY ' 0, 
we have 



V 
with 



y^ti — dct(AT) 



2(K+2 -'ylK+i -r}\f 



y eq2{eJ)JNj 



(12) and eqi{x) 



cos(a:) sin(2:) — a: 



and eq2{x) 



two well 



defined functions on ]0,7r[ that admit well controlled Taylor 
expansion around 0, JiVj = {pt^+2 - Pt ^[a] being the 
skew 3 by 3 matrix (i.e. fcj^^j* — — fc[A]) such that fcj^^j • u = 
fc Am Vfc,u e M.^. 

Special case: rj G Support{T) , ^ T : 



^ eq,{6j)\p,^^,-p,^^,?Nf -Nj 



cos{eJ)eq,{eJ)Nf ■ Nj 

^ 2{\pt,^,-Ti\\pt,^,^^\y 



(15) 



with 6^2(2;) 



,eqi{x) 



cos(a;) sin(a;) — a; 

sin(x)'^ 



p and 

sq3{x) — ''°^^^ly2^ being functions well defined on ]0, 7r[ and 



cos(a:) — 1 

that admit controlable Taylor expansion around 0. 



B. Expression of the Hessians 
We note 5=" = \ ] , (5^ = 




Hwi 



Nf-dy{JmT) 
\ Nf ■d,{Jm'^) 
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det{AT) 



3 



with 



2(K-+2 -v\\pt. 
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■m,,.\.)Nj-Nf.JNj 



3 + 1 

eq^[eJ)d,{Nj)-Nj' -JNJ 



2(bt 



eq,{eJ)Nj .d^jNjf -JNJ 



3egi(^J)(r? - Pt,+2)(c)^7 • Nf ■ JNj 



'ilPt, 



2{\pt^^,-v\\Pt,^,-v\y 

iv - Pt,+-,){c)Nj ■ (2?7 - Pt,+i - Pt.+^y 

K+2 - jjPK+i - ^1^ 

- Ph+2){c)N[ ■ {2r] - pt.^, - pt,^J* 

K+2 - - 

eg8(gJ)(JiVr-Ar7)(,)JiV/ 

2(|Pt,+2-^l|Pt,+x-??l)' 
eqg{ej){2r] - pt,^, - JjVf 

2(K. + 2-^IK + X-^l)' 

(?7-Pt.+x)(c)eg2(^J)JiV7 



2|Pt 



?7p 



{r}-pt,+,)(c)eq2{eJ)JNj 
2^+2 -»?| 



AT/ • 5^* 



K+2 + l 
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and eqeix) = '^^^^^{x)cos{x)/sin{x),eq'^{x) = 
^:^^p-{x)sin{x),eqs{x) = ^:^^^{x)cos{x) / sin{x), and 
eq(){'.i^ = {x)sin{x) being functions well defined on 

]0,7r[ and that admit controllable Taylor expansion around 0. 

Special case: r] G Support{T), ^ T : 



■2\T\^dwf = 



((Pt,+2 - PU+^ Y ■ {Ptj+2 - Ptj+i))(2?7 - Pt,+2 - Pt,+l) 
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Pti + 1 



e9i(^J)((Pt.+2 -ft.+J* • {Pt,^.-Pt,^,))JNj' ■ K 



n\pt,+.-v\\pt,^^~v\r 
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K.+2 -Pu+^lHNf ■ ^n(2r? -pt,+2 -Pt,+J 



2(K+2-'7lK+i-^l)' 



2(b: 



4(bt,+2 -'7lbt,+i 

,71)5 



v\\Pt 



cos{eJ)eq3i9j){JNf ■ Nf + JNf ■ Nj) 



2(K+2 ~v\\Pt_ 



(1 - 2cos(gJ))(iVf ■ Nj){2^-pt,^,-pt^^,) 
2(K.+2-^IK+i-^l)^ 



2 (K+2 -»7lbt 



and e9i(a;) = sin(4'' '^^^(a;) 

2 cos(a:) sin(x)^+3(sin(a;) cos(a:)— x) / \ cos(x) — 1 

sin(x)5 je^s^a;; — sMxy ' 

g„ (2;) = cos(g) sin(a:)''(l-2 cos(x))-2 cos(a)^+2 cos(x) |^gjjj„ 

functions well defined on ]0, 7r[ and that admit controllable 
Taylor expansion formula around 0. 
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(19) 



and 



H{wJ){ri) = ydwJ{n)-n*T 
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